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We present an update on the QCD equation of state of the Wuppertal-Budapest Collaboration, 
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staggered fermion action is utilized. We discuss partial quenching and present preliminary results 
for the fully dynamical charmed equation of state. 
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Figure 1: Continuum limit estimate for the Nf ^2+1 QCD equation of state. No tree level improvement 
was applied to the data points prior to computing the continuum estimate. Perfect agreement is found with 
the Hadron Resonance (HRG) model using physical parameters. Further data is needed for a complete 
continuum limit result. 

1. Introduction 

The properties of the Quark-Gluon- Plasma (QGP) are the focus of a large number of heavy ion 
experiments, such as ALICE at LHC, CERN SPS, RHIC at BNL, and, in the future, FAIR at GSI 
Darmstadt. Hydrodynamical models, based on the observation that the QGP is well approximated 
by an ideal liquid [3, 4, 5], require as an input parameter the equation of state (EOS) of Quantum 
Chromodynamics (QCD). 

Computing the EOS from QCD directly is possible through simulations of the lattice regu- 
larized theory. Lattice QCD. As algorithms and computers improve, such simulations aiming to 
understand the features of the QGP are reaching unprecedented levels of accuracy both for map- 
ping the phase diagram [6, 7, 8, 9, 10] or studying bulk obervables [1, 2, 11, 12, 13, 14, 15, 16, 17]. 
Here, we provide a status report on our efforts to improve upon our result of [2] for the Nf = 2 + \ 
EOS by performing a controlled continuum limit extrapolation and by studying the effects of an 
additional sea charm quark {Nf = 2+1 + 1 EOS). 

2. Continuum estimate for the A^/^ = 2 + 1 EOS 

In order to perform a continuum extrapolation for the EOS, we interpolated the results for any 
given Nt using a standard cubic spline ansatz. Furthermore, we augmented our dataset from [2] 
by adding a set of A'f = 12 points, thereby making a spline interpolation possible. We then varied 
the details of the spline interpolation to account for the statistical uncertainties and, by adding or 
removing spline points, to estimate the systematic effects due to this particular inteipolation pro- 
cedure. The continuum limit was then taken with an 0{a^) fit ansatz. The result for the trace 
anomaly is depicted in Figure 1 . While the data points depicted in Figure 1 have been corrected 
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Figure 2: Finite volume study for our EOS. Shown is the trace anomaly for two different volumes at A^, =8. 
Complete agreement between the two volumes is found, extending our evidence in [2], where the same 
outcome was found for Nt =6. 



for tree level effects, no such correction was applied prior to computing the continuum extrapola- 
tion. As already discussed in [2], these improvement factors would not influence the continuum 
extrapolation procedure. This is also nicely visible in Figure 1 , since the tree level improved data 
points agree perfectly with the continuum estimate. Furthermore, excellent agreement of the lattice 
data with the prediction of the Hadron Resonance Gas (HRG) model was found, strengthening our 
assessment that the results of [2] are accurate. 

Finally, in order to check for finite size effects, we have repeated our analysis at A', = 8 with 
two different aspect ratios, shown in Figure 2. As in [2], where the same analysis was done for 
A', = 6, the results agree perfectly. Therefore, we believe our volumes are more than adequate and 
our results are free of finite size effects. 



3. Partially quenched and dynamical charm effects on the equation of state 

A straightforward approach in estimating charm quark effects is to include the charm quark in 
the valence sector only. Estimates of charm quark effects computed with this ansatz are discussed 
in [2, 19, 20]. Here, the charm quark becomes relevant already for temperatures around 200 MeV, 
which is incompatible with perturbative estimates [18], see Figure 3. 

This is analyzed more closely in Figure 4. Here, the gauge action is plotted in terms of the 
temperature as defined through the LCP for both Nf = 2 + 1 and A'j = 2 + 1 + 1 sea quarks. Note 
that the gauge action enters into the calculation of the equation of state through the integral tech- 
nique (here, p refers to the pressure, and ()*"^ to subtracted/renormahzed quantities; our notation 
is explained in [2]): 

p{T) p{To) .A [^^'"'"'^ ( ,0, , \sub , /,T. .on 
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Figure 3: Left panel: Perturbative calculation [18] of the pressure normalized by the temperature. In this 
perturbative estimate, the three and four flavor results start to deviate for temperatures around 300 Me V. Right 
panel: Summary of different calculations of the deviation ofNf = 2+\ + \ EOS from the Nf = 2+1 EOS, 
both partially quenched and fully dynamical. From top down these results are: [19], [18] (corresponding 
to the result on the left panel), results from a partially quenched calculation including a "Naik" term, a 
continuum estimate for the partially quenched charm contribution, and a fully dynamical calculation at 
Nt = 8, compared to the Nf = 2+1 continuum estimate. Clearly, the partially quenched results consistently 
deviate already at low temperatures, whereas the dynamical results and the perturbative estimate both start 
to deviate at T - 300 MeV. 
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Figure 4: Subtracted negative gauge action ({—Sg) in the language of [2]), used as input in the integral 



technique for calculating the EOS. Shown is the gauge action from the Nf = 2-\ 
partially quenched calculation, compared to the preliminary results for Nf = 2- 
units, no sizable effect is visible. 
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Figure 5: In terms of physical units, the effects of the dynamical charm quark are small. Left panel: 
Subtracted strange condensate. Right panel: Quark number susceptibilities. 

In terms of physical units, no sizable effects from the charm quark in the sea are visible, as is the 
case for other quantities as well (see Figure 5). 

The large deviation between the results shown in Figure 3 is, therefore, likely due to a shift in 
the LCP, which enters into eq. 3.1 through the measure. It is conceivable to try to coiTcct for this 
shift using wNf = 2 + \ + \ LCP, but in the full result some effects will remain. Hence, we prefer 
to perform a fully dynamical simulation of the charmed EOS. 

4. Charmed equation of state for QCD 

Using a new Nf = l + \ + l LCP, we computed the EOS for QCD for A^, = 6, 8, and 10, as 
shown in Figure 6. As already mentioned earlier, the fully dynamical charmed EOS agrees with 
the Nf = 2+ \ EOS up to higher temperatures than the partially quenched one. Similarly to the 
perturbative estimate, the temperature where charm effects become sizeable is close to 300 MeV. 

Furthermore, there do not seem to be sizeable discretization artifacts due to the heavy charm. 
At low temperatures, where the lattice spacing is coarse, the data points for the different Nt agree 
with each other, as well as with the Nf = 2 + \ EOS. As the temperature increases, and the lattice 
spacing becomes increasingly finer, potential discretization effects should become smaller. At our 
present level of precision, we do not see sizable deviations between the different A', from low to 
high temperatures; therefore, we presently believe that discretization effects due to the heavy charm 
quark will not be significant, or at least minor. 

5. Conclusion 

We have presented a first attempt to provide a continuum extrapolated equation of state for 
QCD with Nf = 2+ \ flavors of quarks, and have shown that it is free of finite size effects. We 
intend to supplement this calculation with additional Nt = 12 and very few N = 16 data points. 
In addition, we have addressed the issue of partially quenching for the (charmed) A^y = 2 + 1 + 1 
equation of state, and have argued that the difference of partially quenched and fully dynamical 
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Figure 6: A^/ = 2 + 1 + 1 EOS. The data point are tree-level corrected and compared to the Nf = 2^1 
continuum estimate. The A'/ = 2 + 1 + 1 and Nf — 2~^\ EOS agree well up to a temperature of about 
300 MeV. 

results is caused by a shift in the line of constant physics. Finally, we have shown first fully 
dynamical results for the charmed equation of state. 
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